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ABSTRACT 

We investigate the environments and clustering properties of starburst galaxies se- 
lected from the 2dF Galaxy Redshift Survey (2dFGRS) in order to determine which, 
if any, environmental factors play a role in triggering a starburst. We quantify the 
local environments, clustering properties and luminosity functions of our starburst 
galaxies and compare to random control samples. The starburst galaxies are also 
classified morphologically in terms of their broad Hubble type and evidence of tidal 
merger/interaction signatures. We find the starburst galaxies to be much less clus- 
tered on large (5-15 Mpc) scales compared to the overall 2dFGRS galaxy population. 
In terms of their environments, we find just over half of the starburst galaxies to re- 
side in low to intermediate luminosity groups, and a further ^30 per cent residing in 
the outskirts and infall regions of rich clusters. Their luminosity functions also differ 
significantly from that of the overall 2dFGRS galaxy population, with the sense of 
the difference being critically dependent on the way their star formation rates are 
measured. In terms of pin-pointing what might trigger the starburst, it would appear 
that factors relating to their local environment are most germane. Specifically, we find 
clear evidence that the presence of a near neighbour of comparable luminosity/mass 
within 20kpc is likely to be important in triggering a starburst. We also find that 
a significant fraction (20-30 per cent) of the galaxies in our starburst samples have 
morphologies indicative of either an ongoing or recent tidal interaction and/or merger. 
These findings notwithstanding, there remain a significant portion of starburst galax- 
ies where such local environmental influences are not in any obvious way playing a 
triggering role, leading us to conclude that starbursts can also be internally driven. 

Key words: surveys - galaxies: starburst - galaxies: formation - galaxy: active - 
galaxies: evolution - galaxies: interaction 



1 INTRODUCTION 

Identifying and understanding the drivers of galaxy evo- 
lution remains one of the most important challenges in 
modern astrophysics. Whilst the model of hierarchical large 
scale structure and galaxy formation provides a useful broad 
framework within which to understand how galaxies form 
and evolve, the impact of the environment on the evolution 
of a galaxy remains an important detail that is yet to be 
fully understood. One approach to gauging environmental 
influences on galaxy evolution is to obtain large samples of 
galaxies in the throes of violent evolution enabling searches 
for statistically significant correlations between the galaxy 
properties and the environment it lives in. Thanks to large 
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spectroscopic surveys such as the 2dF Galaxy Redshift Sur- 
vey (2dFGRS; Colless et al. 2001) and the Sloan Digital Sky 
Survey (SDSS; York et al. 2000), such large samples are now 
accessible. Starburst galaxies (galaxies with star formation 
rates orders of magnitude above their quiescent rates) pro- 
vide an excellent opportunity to study galaxies during rapid 
evolution, since the burst is generally expected to have been 
triggered within the last 500 Myrs or so. 

Starburst galaxies are found across the mass spec- 
trum of galaxies, from the massive Ultra Luminous Infrared 
Galaxies (ULIRGs) to dwarf starbursting galaxies such as 
the blue compact dwarfs (BCD) and dwarf irregulars (dlrr). 
There are a number of physical mechanisms proposed for 
triggering star formation across this broad mass spectrum, 
some of which are due to intrinsic processes, whilst others 
are due to environmental interaction. Spiral density waves 
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are known to trigger star formation, however this mechanism 
is probably not strong enough to trigger a massive starburst 
in larger galaxies, whilst bars are effective at transporting 
gas to the ce ntral regions of a galaxy triggering a circum- 
nuclear burst l|Kennicuttl 19981 ) . Whilst inadequate or absent 
in massive galaxies, internal pr ocesses such as stochastic self- 
prop agating star formation llGerola, Seiden fc Schulmannl 
|198C| ) and cyclic gas re-processing jpapaderos et al.lll996l ) 
may cause dwarf galaxies to undergo multiple cycles of star- 
bursts and quiescent phases on timescales that are short 
compared to a Hubble time. 

Major mergers of two or more gas-rich galaxies have 
been shown to produc e massive bursts of star forma- 
tion (|Bekki et al.l l2001al ; iMihos fc Hernquistl Il996l ), whilst 
more than 95 per cent of ultra-lu minous infrared galaxies 
(ULIRGs) are merger systems (see lSander fc Mirabell (|l996l ) 
for a review). Numerous studies have been conducted on 
the effects of close galaxy pairs on star formation rate 
(SFR), with the general consensus being pairs with pro- 
jected separations of less than ~ 5 kpc show enhanced 
SFRs ilFreedman Woods et all 120061 ; iNikolic et ail |2004 
lLambas et al.ll2003h . There is evidence that the mass ratio 
of the pairs may also be important. iFreedman Woods et al] 
(2006) find evidence that pairs with |Amn| < 2 show strong 
correlations between SFR and projected separation, whilst 
those with |Am.R| > 2 show no correlation, consistent with 
lLambas et ail (|2003h . These results suggest that tidal inter- 
actions and mergers with galaxies of approximately equal 
mass are important in triggering star formation. 

The effects of the cluster environment may trigger 
starbursts; theoretical studies have shown that it is possi- 
ble to trigger a starburst due to the high static thermal 
and ram pressure felt by molecular clouds within a galaxy 
disc due to the tenuo us, hot intra-cluster medium (ICM) 
(Bckk i~fe Couchll2003l ). whilst group member galaxies may 
experience starbursts during the infall process due to th e 
time dependant tidal gravitational field (|Bekki et al.ll 19991 ). 
The above mechanisms can be tested observationally, since 
the different triggering mechanisms are expected to produce 
different spatial distributions of the starburst across the 
galaxies. Much observational work has been carried out on 
the fractions of the types of galaxies residing in clusters, and 
on the populations in major cluster mergers. Photometric 
studies revealed significant evolution of the blue population 
with r edshift (the Butcher- Oemler effect . iButcher fc Oemlen 
(1978)). Subsequent spectroscopic studies revealed the effect 
was mostly caused by galaxies which had undergone a burst 
of star formatio n between 0.1 and 1.5 Gy rs prior to the epoch 
of observation (|Couch fc Sharpleslll987i ). 

At the low mass end of starburst galaxies are the BCDs 
where much work has be en carried out in order to determine 
triggering mechanisms. iPustilnik etUI (|200ll ) find that a 
significant fraction (~ 80 per cent) of a sample of 86 BCDs 
either have neighbours close enough to induce strong tidal 
forces, or display a disturbed merger morphology, imply- 
ing that i nteractions are impor tant in triggering starbursts 
in BCDs. iNoeske et alj (1200 ll ) draw a similar conclusion 
from a sample of 98 star forming dwarf galaxies, finding a 
lower l i mit of ~ 30 per cent harbour close neighbours, whilst 
iTavlorl (|l997h find the rate of companion occurrence in HII 
galaxies is ~ 57 per cent, more than twice the rate of com- 
panions in a comparison sample of low surface brightness 



galaxies. However, iBroschl {2004) find that in a sample of 
96 BCG and dlrrs imaged in Ha there is not a significant 
number of star forming companions, concluding that tidal 
triggering is probably not important in enhancing star for- 
mation in these galaxies. Based on the observation that the 
star formati on seems to be occurring at the outer boundaries 
of the dlrrs, IBroschl (|2004l ) propose a 'sloshing' mechanism 
for triggering the star formation, where gas oscillates within 
the dark matter halo, producing gas build up at locations 
where the gas 'turns around' providing a trigger for star 
formation. In any case, the triggering mechanism for these 
dwarfs is still under much debate. 

Also of relevance in this context is the population of 
post-starburst galaxies, first discovered in distant clusters 
and called 'E+A' galaxies by Dressier & Gunn (1982, 1983). 
These galaxies are known as post-starburst galaxies, since 
their rather enigmatic spectra - characterized by a strong A- 
type stellar spectrum superimposed upon a normal E galaxy 
type spectrum and devoid of any emission lines - would in- 
dicate they have undergone a recent burst of star formation 
that was abruptly truncated within the last 1 Gyr (Couch 
& Sharpies 1987, Poggianti et al. 1999). Subsequent stud- 
ies have shown that E+A galaxies are also found in a va- 
riety of environments at lower redshifts (e.g. Zabludoff et 
al. 1996), and o f parti cular note is the 2dFGRS-based study 
by iBlake et~all (|2004T l. They identified all the E+A galax- 
ies within the 2dFGRS and studied their local and global 
clustering environments in an attempt to pin-point what 
might have triggered the original starburst. They concluded 
that the environments of E+A galaxies do not differ signif- 
icantly from that of the general 2dFGRS galaxy population 
as a whole. However, the morphologies are consistent with 
a merger origin, and hence the triggering/cessation mecha- 
nism for the initial burst is likely to be driven by very local 
encounters. Simulations have shown that the E +A signature 
marks the late stage of a galaxy-galaxy merger (|Bekki et al.l 
l2001bl ) where the merging companions are no longer sepa- 
rable. The next logical step, therefore, is to study the envi- 
ronments of the starburst galaxies, which are thought to be 
the progenitors of E+A galaxies. 

Motivated by this need and the broader issues outlined 
above as to what mechanisms are responsible for triggering 
starburst events, this paper reports a systematic study of 
the environments and clustering properties of the starburst 
galaxies within the 2dFGRS, tha t very much parallels the 
E+A study bv lBlake et al.l (120041 ). We select two starburst 
samples, using complementary star formation rate determi- 
nants, with each sample containing 418 galaxies. This al- 
lows measurement of statistical properties of the starburst 
population, such as the luminosity function, near neighbour 
properties, morphological distribution and large scale clus- 
tering. The paper is organised as follows: In Section [2] we 
briefly introduce the 2dFGRS and describe the selection 
criteria for our starburst samples. In Section [3] we use the 
Supercosmos Sky Survey (hereafter SSS) images to study 
the morphologies of our samples. In Section [4] we measure 
the 6j band luminosity function. In Section [S] we inves- 
tigate the small scale environment surrounding the star- 
bursts. In Section [6] we measure large scale environmental 
properties of the sample. In Section [7] we summarise and 
discuss our findings. Throughout this paper, we assume a 
standard ACDM cosmology where Sl m = 0.3, Q,a = 0.7 and 
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h = tfoAlOOkms" 1 Mpc" 1 ) = 0.7. We convert redshifts to 
physical distances using these parameters. 



continuum level from effects such as scattered light. We cor- 
rect for broadening of the spectral lines due to redshift by 
dividing the EWs by 1+z, where z is the galaxy redshift. 



2 SAMPLE SELECTION 

2.1 The 2dF Galaxy Redshift Survey 

Our starburst samples were selected from the final release 
of the 2dFGRS, a large spectroscopic survey conducted at 
the Anglo- Australian Telescope using the Two-Degree Field 
(2dF) multi-fibre spectrograph. The 2dF instrument is capa- 
ble of simultaneously observing 400 objects over a 2° diam- 
eter field of view. The survey consists of ~220,000 galaxies 
covering an area of approximately 2,000 deg 2 in three re- 
gions: a South Galactic Pole (SGP) strip, a North Galactic 
Pole (NGP) strip and a series of 99 random fields dispersed 
around the SGP. The input source catalogue for the 2dF- 
GRS was derived from an exte nded version of the Auto- 
mated Plate Measuring (APM; IMaddox et al.lll990a. b. d ; 



iMaddox. Efstathiou fc Sutherland! 1996T ) galaxy catalogue 
based on APM machine scans of 390 plates from the UK 
Schmidt Telescope (UKST) Southern Sky Survey. The ex- 
tinction corrected survey limit was &j = 19.45 and the me- 
dian redshift is z = 0.11. The spectra are collected through 
2 arcsec fibres, cover the wavelength range 3600 - 8000A, 
have two-pixel resolution of 9.0 A and a median S/N of 
13 pixel" 1 . The large range in wavelength is made possi- 
ble by the use of an Atmospheric Dispersion Compensator 
(ADC) within the 2dF instrument. Redshift determinations 
are checked visually and allocated a quality parameter, Q, 
ranging from 1 to 5. Redshifts with Q ^ 3 are regarded 
as reliable at the 98.4 per cent level. The rms uncertainty 
of the redshifts is 85 km s _1 . The spectra are not flux cal- 
ibrated and consist of a sequence of 'counts' as a function 
of wavelength. The 2dFGRS data base can be accessed at 
http:// www.mso.anu.edu.au/2dFGRS 



2.2 The 2dFGRS spectral line catalogue 

We used the 2dFGRS spectral line catalogue prepared by 
Ian Lewis to select o ur samples. The fit ting procedure is 
described in detail by iLewis et all (2002). Here we briefly 
outline the fitting procedure as follows: The continuum in 
the vicinity of each spectral line was removed by subtract- 
ing the median level over a 133 A window. Gaussian profiles 
were fitted for up to 20 lines per spectrum in both emission 
and absorption, with the height, width and a small pertur- 
bation in redshift (always Az < 0.005) as free parameters. 
After fitting, the quality of the fit was determined using the 
parameters obtained from the fit (height, width and area) 
and the rms residuals. A quality flag was assigned to each 
line fit in the range (bad fit) to 5 (good fit) and a signal to 
noise parameter was calculated for each line by taking the 
average signal per pixel and dividing by the average noise per 
pixel for pixels lying within ±3cr of the line center. Equiva- 
lent widths (EWs) were calculated using the measured flux 
in the line and the fitted continuum level (in the 133 A win- 
dow local to the line). The EWs can be used in line ratio 
diagrams without the need for flux calibration. However, it 
is assumed there is no significant flux contribution to the 



2.3 Selection of emission line galaxies 

In order to select a high quality sample of starburst galax- 
ies, we first ensured the parent catalogue consisted of single 
observations of galaxies with high quality spectra. Thus, we 
kept only spectra with: 

• ADC=lQ 

• 2dFGRS redshift parameter, Q ^ 3 

• z ^ 0.002 to exclude stars. 

The parent 2dFGRS spectral line catalogue contains 
264,765 spectra, and the exclusion of those objects failing 
the above criteria reduced this to 172,427 spectra. We then 
included only the highest quality spectra of the multiply ob- 
served objects based on the redshift quality flag or the high- 
est mean S/N measured between 4000 and 7500 A (if the 
redshift quality flags were equal) , leaving a parent catalogue 
of 162,529 high quality spectra with no multiple observa- 
tions. 

From this catalogue, we selected a sample of emis- 
sion line galaxies from which a sub-sample of star form- 
ing galaxies was selected. It is standard practice to dif- 
ferentiate star forming galaxies from narrow-emission-line 
AGNs according to their position on the BPT diagram 
^Baldwin. Phillips fc Terlevicbl [l98lh . where the logarithm 
of the ratio's of [OIIIJ/H/3 and [Nil] /Ha are plotted. With 
this in mind, we selected galaxies with EWJS 2.5 A in emis- 
sion, line S/N> 1.0 and a line fit quality flag ^ 4 for the fol- 
lowing species: Ha, Nil, OIII and H/3. The cut of EW^ 2.5 A 
is arbitrary, but ensures a robust line detection for use in the 
line ratio diagram. We detected 9,516 emission line galax- 
ies that conform to the above criteria; Figure 1 shows them 
plotted in the BPT diagram. The advantage of the chosen 
lines for the ratios plotted is three-fold. Firstly, the lines are 
close together in wavelength space, hence extinction due to 
dust has little effect on the calculated ratios. Secondly, the 
change in the continuum level is minimal for nearby lines, 
meaning the ratio is essentially equivalent to a line flux ratio. 
Thirdly, as seen in Figure 1 the plot of the ratios show two 
distinct populations in the [OIII]/H/3- [Nil] /Ha plane, the 
origin of which lies in the different physical mechanisms pro- 
viding the photo-ionizing continuum. For AGNs, the ionizing 
radiation is non-thermal in nature and is well approximated 
by a power law, whilst the ionizing radiation in star-forming 
galaxies is in t he form of UV radiatio n emit ted by hot OB 
type stars (see IVeilleux fc Osterbrockl (|l987T l for a detailed 
description). The empirically derived starburst/ AGN classi- 
fication line of Kauffmann et al. (2003) was used to select a 
sample of galaxies which have emission dominated by star- 
formation (although we note there may be a small contri- 
bution to emission line strength due to AGN emission) such 
that a star forming galaxy is one which has 



1 Prior to August 1999, the ADC suffered positioning problems 
resulting in the blue end of each spectrum being severely depleted 
of counts. All spectra obtained during this period have been as- 
signed ADC flag =0 
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Figure 1. BPT emission line ratio diagnostic diagram, with all 
the 2dFGRS emission line galaxies plotted as dots. The solid line 
represents the empirically derived star formation/ AGN demarca- 
tion line. Due to the non-thermal mechanism for the underlying 
photo-ionising continuum, AGNs inhabit a different region of this 
diagram compared to star-forming galaxies, as can be seen by the 
plume of galaxies to the upper right. Galaxies residing in this 
plume are excised from our sample as their emission is not due to 
star formation. 
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This sample contains 8,369 galaxies from which our star- 
burst samples are selected. 

It is pertinent to note here that the 2dFGRS fibres 
cover only a 2-2.16 arcsec angular diameter on the galaxies 
(4kpc at the median 2dFGRS redshift of z = 0.11), and 
our sample will not contain galaxies where star formation is 
not occuring within the area covered by the fibre aperture. 
However, the mean RMS error in the fibre positioning was 
0.3 arcsec, whilst the seeing was of the order 1.5-1.8 arcsec, 
hence the observed spectra cover a larger fraction of the 
galaxy once these effects are accounted for, and the aper- 
ture effect is somewhat diluted. We also note that since we 
do not obtain the total Ha flux of each gal axy, accurate 
total star formation rates cannot be derived. iKewlev et al.l 
(2005) showed that properties such as total aperture cor- 
rected SFRs derived from 2dFGRS spectra wou ld be most 
affected at z < 0.06, whilst [Hop kins et al.l (|2003l ) found that 
those galaxies with the largest aperture corrections (ie. those 
with the largest angular size) had their Ha SFRs slightly 
overestimated when compared to the 1.4 GHz SFRs. Thus, 
we expect galaxies at lower redshifts will have their aperture 
corrected Ha SFRs slightly overestimated. In section 12.4.21 
we attempt to correct for this and other redshift dependen- 
cies by binning the data in redshift space when selecting our 
starburst sample. The combination of redshift binning and 
the selection of only the most extreme galaxies in SFR per 
redshift bin should serve to minimize the aperture effect, 
although we note there may be some contamination in our 
sample. 

Figure [2] shows that the Ha EW (hereafter Who) does 
not show any trends with redshift in the range 0.04 < z < 
0.16. However, we find evidence of a slight increase in the 
lowest redshift bin (0.002 < z < 0.04) with the median 
Whu increasing from 40 to 60. This increase is due to 
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Figure 2. Plot of Wua as a function of redshift. The solid lines 
show (from bottom to top) the first, second (median) and third 
quartiles. They show that there is no appreciable trend with red- 
shift, apart from at low redshifts where the turn-up in Wjja is 
possibly due to contamination by extra-galactic HII regions and 
dwarf starbursting galaxies that are not detected at higher red- 
shifts. The gaps in the distribution are due to redshifting of emis- 
sion lines used in the selection process into areas where sky dom- 
inates the spectrum. This causes the lines to have bad fits, and 
thus be rejected by our criteria. 



low mass emission line systems such as HII regions and 
dwarf starbursting galaxies which, due to their low luminos- 
ity, would not be detec ted at higher redshifts. Interestingly, 
iMateus fc Sodrd l|2004f ) find no trends in Whc with redshift 
out to z=0.05. The gaps in the Who distribution correspond 
to redshifts where emission lines used in the selection crite- 
ria are shifted into areas where sky emission/absorption lines 
dominate the spectra, producing bad line fits meaning these 
galaxies are rejected during the selection process. Whc is a 
relative quantity meaning the effect of aperture depends on 
both the spatial distribution of the Ha line strength and the 
continuum profile in the galaxy. We assume this distribution 
is uniform and we make no attempt to correct for aperture 
in the starburst sample in Section ^. 4. II We note there will 
be contamination of our sample due to HII regions, however 
this contamination is negligible and makes up only 1.6 per 
cent of the sample selected in Section 12.4.11 



2.4 Starburst galaxy selection criteria 

Starburst galaxies are generally defined by their extremely 
high SFRs, such that if the SFR was sustained, the gas 
reservoirs feeding the burst will be exhausted in timescales 
much shorter than the Hubble time, typically of the order 
of 10 s yrs. The problem then becomes selecting observable 
quantities that adequately fit the above definition. In gen- 
eral, an observable quantity is converted into an SFR via 
some theoretical model which makes various assumptions. 
Some examples of these observables are: Ha luminosity, in- 
frared luminosity, radio emission, a galaxies position on the 
[O H1-H<5 plane, [OH] luminosity, UV lum i nosity etc (see 



eg. ISullivan et al.l |200ll ); iKennicuttl |l99Sl) ; iHopkins et all 

(2003)). A galaxy is then generally defined as a starburst 
based on some arbitrary cutoff in absolute SFR, SFR per 
unit area or the Scalo b-parameter (ie. the ratio of current 



Environments of 2dFGRS starbursts 5 



to average past SFR). Below we describe two samples se- 
lected from the emission line galaxies above using one ap- 
proach equivalent to measuring the b-parameter and one 
equivalent to measuring the absolute SFR. 



2.4-1 Normalized SFR sample 

The first method employed to estimate the SFR of the 
emission lin e gala xies selected above uses the method of 
I Lewis et alJ l)2002h to determine a normalized SFR, (i* , de- 
fined by 

H*=r,W Ha L*, (2) 

where r\ = 8.2 x f0~ 4 2 Mpsyr~ 1 erg~ 1 i s the conver- 
sion from L Ha to SFR (|Sullivan et al.ll200ll ), L* = f.f x 
10 40 ergs _1 A 1 is the luminosity corresponding to the 
knee in the r' band luminosity function as determined by 
iBlanton et all \2004 ) and Whc is the equivalent width of 
the Ha emission line (Whc ~ Lhc/ L c , L c is the continuum 
luminosity) . The SFR measured in this way gives a SFR nor- 
malized to L* . A 1 A c orrection is made to Whc to account 
for stellar absorption l|Hopkins et alj|2003h . Thus the SFR 
normalized to L* is given by fi* — 0.0902(W / £f a + f). 

It is not clear where the dividing line between a nor- 
mal star forming galaxy and a starbursting galaxy can be 
drawn using this method for determining the SFR, hence 
we arbitrarily defined our starburst sample to be the top 5 
per cent of galaxies in terms of their /«* values; hereafter we 
refer to this as our SFRNORM sample. This corresponds to 
a cut at (i* = 9.6M Q yr _1 , or a W Ha of ~ 106 A. W Ha is 
essentially a measure of the ratio of the SFR to the mass of 
the old stellar popul ation, and it correlates strongly with the 
scalo b-parameter (|Kennicutt et al.1ll994l ). The lower limit 
on b correspond i ng to our Who is ~ 2 (using Table 1 of 
iKennicutt et alJ (|l994T ) with a Salpeter IMF). Since qui- 
escent disks have b ~ 0.5, th is seems a reasonable cutoff 
value to delineate starbursts (|Kennicutt et al.1 [2005). Fig- 
ure [3] shows example spectra selected from the SFRNORM 
sample where the spectra are de-redshifted and normalized 
such that the strongest emission line peaks at 1. Dominant 
emission species are marked on the plots. 



where (Who) is corrected by lA to allow for the neg ative 
contribution of stellar absorption (|Hopk ins et all [2003) and 
the factor of 2.5 corrects for ~ 1 magnitude extinction of 
the tab magnitude due to dust exti nction, con s istent with 
the extinction measured at Ha b y iKennicuttl (Il983l ) and 
iNiklas. Klein fc Wielebinskl (|l997t ). The star formation rate 
(SFR) is related to L Ha by eqn. 1 of ISullivan et al.l i|200ll ); 

SFR(M e yr" 1 ) = 8.2 x HT 35 L Hct (W). (5) 

The discrepancy between the effective wavelength of the tab 
band and the wavelength of the Ha emission should not have 
a large effect on our measured SFRs since the gradient of the 
continuum at these wavelengths is close to zero. This method 
effectively corrects for aperture effe cts, assuming the Wh c, 
is constant across the galaxy (see iHopkins et al. 

I <|2003l )). 

hence it is equivalent to measuring the total absolute SFR, 
hereafter referred to as SFRTOT. 

Figure [4] shows the dependence of the SFRTOT values 
on redshift. There exists a strong correlation with redshift, 
part of which can be explained by a combination of the sur- 
vey magnitude limit, and a lower cut-off in Whc at around 
20 A. This produces an effective limit in the SFR observed 
with redshift. It is also known that below z ~ 0.06, aper- 
ture effe cts cause discrepan t measures of SFR u sing this 
method l|Kewlev et al.1 120051 ; IHopkins et~ai1 120031 ). We at- 
tempt to nullify any redshift dependency by binning the 
sample in redshift space, such that each redshift bin con- 
tains 500 galaxies (to ensure equal statistical weight), and 
selecting the galaxies in the top 5 per cent of the SFRTOT 
distribution within each bin. This serves to select the most 
extreme star-forming galaxies in each redshift bin, and hence 
defines our second ('SFRTOT') starburst sample. Figure [4] 
shows the cutoff in SFRTOT for each bin, where again the 
gaps correspond to regions where emission lines are red- 
shifted into sky dominated portions of the spectrum, hence 
are rejected by our selection criteria due to bad line fits. 
Example spectra of a sub-sample of galaxies selected from 
the SFRTOT sample are presented in Figure [5] in a similar 
fashion to Figure [3] As for the SFRNORM sample, these 
spectra are clearly dominated by emission lines produced 
both by photo-ionization and forbidden transitions. 



24.2 Total SFR sample 

It is common to flux calibrate spectra using broad band pho- 
tometry. We used tf magnitudes obtained from the SSS to 
estimate the continuum level at the wavelength of Ha. The 
tf magnitudes were K-corrected to th eir rest frame values 
using equation 37 of I Wild et all (|2004l ) and the K-corrected 
values were con verted to tab ma gnitudes with equations 
A12 and A16 of ICross et al.1 l|2004l ). We use equation 15 of 
IBlanton et aD (|2004l ) to derive the continuum luminosity, 



: 4.18 x 10 24 (WA _1 ) x 10 



-0.4M A i 



(3) 



where Mab is the absolute AB magnitude and the corre- 
sponding effective wavelength for the tab band is A e // = 
5595 A. The Ha luminosity is then derived from L cont and 
the Ha equivalent width such that 



2.5L cont x (W Ha + 1), 



(4) 



2.4.3 Pristine SFRTOT starburst sample 

The above samples select out starburst galaxies based on the 
SFR derived from Ha alone. The advantage here is that Ha 
is less affected by dust than bluer lines such as H<5 and [Oil] 
- both common tracers of recent star formation. Thus in 
the above, there is no restriction placed on the H<5 and [Oil] 
lines. This means we are not biased against selecting 'dusty' 
starbursts which have their emission at bluer wavelengths 
heavily attenuated. However, this also means we may be 
selecting out older starburst galaxies which have been form- 
ing stars at a high rate for a several hundred million years. 
Indeed, it was noticed that a high proportion of the SFR- 
TOT sample had H<5 in absorption. This may mean we are 
selecting a combination of both dusty starbursts and older 
starbursting systems. To test for any effect this may have, we 
have selected a sub-sample from the SFRTOT sample with 
the criteria that H5 is either in emission ( EW^ and qual- 
ity flag ^ 4) or is not detected in absorption or emission (its 
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Figure 3. Example spectra for starburst galaxies within our SFRNORM sample (see Fig. [7] for corresponding SSS images). The 
wavelength scale is in the galaxy rest frame. Prominent night-sky emission lines or absorption features have been excised and interpolated 
over. Dashed vertical Lines show the positions of prominent spectral features. 



quality flag is 0), and hence the nebular emission is domi- 
nating the stellar absorption. This ensures we are selecting a 
sample where the majority of the light e mitted comes from a 
starburst that started ~ 200 Myrs ago (jBalogh et al.ll 19991 ). 
hence we are more likely to find evidence of a triggering 
mechanism. This added selection criteria produces a sub- 
sample of 244 galaxies within the SFRTOT sample, which 
we will call the SFRTOT+ sample. In the following analy- 
sis, we present results for the extra sub-sample in addition to 
the main SFRTOT sample. A sub-sample of the SFRNORM 
sample using the above criteria will not be presented, since 
97 per cent of the SFRNORM galaxies meet the above crite- 



3 MORPHOLOGIES 

Of relevance to our study is the morphology of our differ- 
ent samples of starburst galaxies and, in particular, whether 
there is direct evidence that mergers and tidal interactions 
are responsible for driving the starburst activity. All the 
galaxies in our SFRNORM and SFRTOT starburst sam- 
ples were therefore morphologically classified by one of us 
(W.J.C.) by visual inspection of blue and red (tf) im- 
ages obtained from the SSS. The SSS has digitised sky sur- 



vey plates taken with the UK Schmidt telescope, with a pixel 
size of 0.67 arcsec. The data were accessed using the website 
http:/ /www-wfau. roe. ac.uk/sss/pixle. html A description of 
the SSS is given by Hambly et al. (2001). 

The scheme used to morphologically classify our star- 
burst samples was as follows: The first step was to try and 
assign a basic Hubble type to each galaxy. In doing so, the 
following five broad categories were used: E or SO (E/S0), 
early-type spiral (SE), medium-type spiral (SM), late-type 
spiral (SL) and irregular (Irr). Using these broader cate- 
gories, rather than attempting to use all the individual Hub- 
ble types (Sa, Sab, Sb, etc), was considered to be sufficient 
for our purposes and really all that was practical given the 
limited quality and resolution of the SSS images. The num- 
bers of galaxies classified into each of these five morphologi- 
cal bins are plotted in histogram form in Figure [6] and listed 
as percentages in Table 1, for both the SFRNORM and SFR- 
TOT samples. We see that both samples have a broad mix 
of morphologies, with no one type being dominant. Only in 
the SFRTOT sample do we see the absence of a particu- 
lar morphological class, which somewhat surprisingly is the 
Irr types. Figure [7] shows example fej band SSS images of 
5 galaxies each from the SFRNORM and SFRTOT sam- 
ples (corresponding spectra are presented in Figures [3] and 
[5]). The plots are labelled with the morphological classifi- 
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Figure 4. Plot of the total star formation rates (SFRTOT) measured for the 2dFGRS galaxies versus redshift. The green line indicates 
the effective SFR limit derived from the apparent rp band magnitude limit and WHa ~ 20 corresponding to approximately the lowest 
measured Wjja m the star forming sample. The blue horizontal bars indicate our binning of the data in redshift (to ensure equal statistical 
weight to each bin) with their vertical height indicating the top 5 per cent cutoff for defining our SFRTOT starburst sample. 



cation, absolute bj magnitude, redshift and 2dFGRS label. 
No morphological selection has been applied in selecting the 
starbursts appearing in this figure. 

Whilst the majority of galaxies in both samples could 
be assigned a Hubble type, there were nonetheless a sub- 
set of galaxies where this was not possible. Such galaxies 
fell into one of two cases: (i) They were sufficiently well re- 
solved to discern their structure, but that structure was so 
peculiar that it did not resemble any of the Hubble types. 
These cases are denoted 'Pec' in Figure[S]and Table 1, where 
it can be seen that they represented only 2 — 3 per cent of 
each sample, (ii) They were too small to be resolved and 
hence be classified. Notably, we only encountered galaxies 
of this type - which we denote as 'Comp' (for compact) - in 
the SFRNORM sample, where they account for 16 per cent 
of these types of starburst galaxies. The 'Comp' classified 
galaxies can be separated into two groups: (i) The low red- 
shift (z < 0.05) galaxies whose compact nature, conspicu- 
ously blue colours (as seen on the SSS images), and their 
generally faint, 'dwarf-like luminosities (Mj,j > —18), sug- 
gest they are most likely blue compact dwarf (BCD) galaxies 
jThuan fc Martini fl98ll) . (ii) Those galaxies with z > 0.05 
and M(,j < — 18 where the SSS images are simply not of 
high enough quality to assign a Hubble type. 



The second step of our morphological classification pro- 
cess was to assess the images for evidence of on-going or re- 
cent merger and/or tidal interaction activity. To be classified 
as an on-going merger, a galaxy had to exhibit clear visual 
evidence of it coalescing with another, as manifested by the 
presence of tidal bridges and/or arms. Similarly, to be classi- 
fied as being involved in a tidal interaction, a galaxy had to 
show clear signs of a tidal link (via a bridge or arm) with an- 
other galaxy. As such, we mainly identified what would be 
called 'major' mergers/interactions (where the two galax- 
ies involved have similar luminosities), given that it is only 
these types that are sufficiently conspicuous to satisfy the 
rather conservative visibility criterion that we took. Finally, 
evidence of a recent merger or interaction was inferred from 
the presence of tidal tails, shells or debris. This was by far 
the most difficult and subjective classification to make. In- 
deed, sometimes the evidence for a merger or interaction, 
while being highly suggestive, was not always totally con- 
vincing, and such cases were placed into a second category 
called 'Possible Merger/Interaction'. 

The numbers of starburst galaxies classified as merg- 
ers/interactions in this way are shown by the dark shaded 
histograms in Figure [6] and listed as percentages in the last 
two columns of Table 1. Importantly, we find that a sig- 
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Figure 5. Example spectra of starburst galaxies in our SFRTOT sample, displayed in the same manner as Fig. [3] 
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Figure 6. Histograms of the number of galaxies per morphological type within the SFRNORM (left panel) and SFRTOT (right panel) 
samples. The filled portions of the histograms show the number of galaxies within that particular class which showed evidence of a merger 
or some form of tidal or interaction activity. 



nificant fraction of our starburst galaxies are involved in 
on-going or recent mergers and interactions. The highest 
number is found in the SFRTOT sample where 28 per cent 
of the galaxies are classified as mergers/interactions, with a 
further 6 per cent possibly being in this category. The same 



percentages for the SFRNORM sample are 16 per cent and 
6 per cent, respectively. While these percentages are signif- 
icant and hence identify the mechanism likely for trigger- 
ing the starburst in somewhere between a fifth and a third 
of the cases, the fact nonetheless remains that the major- 
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Table 1. Morphological classifications for the SFRNORM and SFRTOT samples, based on visual inspection of SSS fej and r-p images. 
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Figure 7. SSS bj band lxl arcmin images of the starburst galaxies whose spectra are shown in Figs. 3 and 5. The top row contains 
the images for the SFRNORM galaxies shown in Fig. 3 while the bottom row contains the images for the SFRTOT galaxies shown in 
Fig. 5. The bar on the lower right of each image corresponds to 10 kpc at the galaxy redshift. The circle shown at the top right hand 
corner indicates the size of the 2arcsec 2dF fibre. The absolute 6j magnitudes are given in the bottom left corner (as Mj, — 51og 10 /i). 
The morphological classifications from Section \3\ are given in the top left corner 



ity (70 — 80 per cent) of the galaxies must have had their 
starbursts triggered in some other way. Here our visual in- 
spection of all our objects provided a possibly important 
additional clue. While many of our starburst galaxies were 
isolated galaxies (no near neighbour within the 1 arcmin on a 
side SSS images used for our morphological classifications), 
the number that were in crowded fields with one or more 
very close neighbours was even more striking. Hence an im- 
portant part of our analysis is to examine the 'local' envi- 
ronments of our starburst galaxies and rigorously quantify 
whether such visual impressions are significant and deter- 
mine whether the presence of a bright or faint neighbour 
is of any relevance in searching for other starburst trigger 
mechanisms (see Section 5). 



4 LUMINOSITY FUNCTION 

The luminosity function is an essential diagnostic tool for 
samples of galaxies. It encodes information about galaxy 
populations which may hold key clues to their evolution. 
For example, if it is found that our population comprises 
only very luminous galaxies we could speculate that in order 
to trigger a starburst a major merger would be required, 
whereas a sample of low luminosity galaxies are more likely 
to be affected by tidal interactions, hence may be triggered 
by weaker interations. 



We derive the foj band luminosity functions, $(M), 
for our three (SFRNORM, SFRTOT and SFRTOT+) star- 
burst samples and the 2dFGRS sample using the Step- 
wise Maximum Likelih ood (SWML) method described by 
lEfstathiou et al.1 (1 19881 ). In order to minimize the effect of 
varying completeness in the different 2dFGRS fields, we set 
apparent magnitude limits of 14 < bj< 19.2 for all 3 sam- 
ples. For the complete 2dFGRS sample, the redshift was 
limited to 0.002 < z < 0.18 to match the redshift lim- 
its of our starburst samples. These limits restricted the 
initial samples to sizes of 311, 412, 239 and 108638 for 
the SFRNORM, SFRTOT, SFRTOT+ and 2dFGRS cat- 
alogues, respectively. SSS 6j and r-p magnitudes were used 
to determine absolute magnitudes which were 'K-corrected' 
and also 'E-corrected' f or the luminosity evo lution of an av- 
erage 2dFGRS galaxy (|Norberg et al.ll2002D . 

The SWML estimator requires normalization, hence we 
normalize each luminosity function to a constant source sur- 
fac e density, a. To obtain a, we used the parameters derived 
by iNorberg et all |2002) in fitting a Schechter function to 
the 2dFGRS catalogue in the absolute magnitude interval 
-16.5 < M bj -51og 10 h < -22 (i.e., -51og 10 h = -19.66, 
a = -1.21 and <F = 1.61 x 10" 2 /i 3 Mpc -3 ). We integrate 
the Schechter function using these parameters within the 
redshift interval 0.002 ^ z ^ 0.18 and magnitude interval 
14 < 6j < 19.2 to obtain a = 124.6 deg" 2 . This normal- 
ization allows us to investigate differences in shapes of the 
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luminosity functions for each sample. The fractional errors 
presented for <I>(M) are Poisson l/\/~N errors, where N is the 
number of objects in the magnitude bin of interest. We used 
the least-squares method to fit Schechter functions to each of 
the 2dFGRS, SFRNORM, SFRTOT and SFRTOT+ lumi- 
nosity functions, with the best fitting Schechter parameters 
presented in Table [2] Figure [8] shows the luminosity func- 
tions derived from the SFRNORM, SFRTOT, SFRTOT+ 
and 2dFGRS samples with best fitting Schechter functions 
overplotted. 



Table 2. Schechter function parameters derived from a least- 
squares fit to the 2dFGRS, SFRNORM, SFRTOT and SFRTOT+ 
samples. The errors are ±lcr. 



Sample 


M b, ~ 51ogi (h) 


a 


(10" 3 h 3 Mpc- 3 ) 


2dFGRS 


-19.65 ±0.02 


-1.17 ± 0.01 


13.1 ± 0.3 


SFRNORM 


-18.73 ±0.20 


-1.50 ±0.07 


27.2 ± 7.6 


SFRTOT 


-18.51 ± 0.07 


1.24 ±0.14 


13.8 ± 1.4 


SFRTOT+ 


-18.69 ± 0.1 


0.58 ±0.16 


21.1 ± 1.4 



The appearance of the SFRNORM luminosity function 
is significantly different from that of the 2dFGRS sample, 
and we quantify this using a chi-squared test. The chi- 
squared test is performed only where bins contain more 
than five objects (since the chi-squared test is not ro- 
bust where bins contain less than 5 counts), and returns 
a vanishingly small probability that the SFRNORM lumi- 
nosity function is consistent with that of the whole 2dF- 
GRS catalogue. This also becomes clear on comparison of 
the SFRNORM and 2dFGRS Schechter function parame- 
ters, with the SFRNORM sample showing both a fainter 
Mj — 51og 10 (h) and steeper faint end slope (a). It is clear 
that the SFRNORM sample is dominated by galaxies with 
Mb 3 — 51og 10 (h) > —18 when compared to the overall 2dF- 
GRS sample. This portion of the luminosity function is gen- 
erally occupied by dwarf and irregular type galaxies. The 
domination of low luminosity galaxies in the SFRNORM 
sample is not unexpected, since the normalized SFR is es- 
sentially a SFR per unit luminosity, hence the measurement 
is biased towards low luminosity galaxies. This bias does 
not detract from the physical importance of this sample; 
the SFRNORM sample just consists of a different class of 
starburst object, i.e., those that don't necessarily have prodi- 
gious SFR, but have SFR much higher then their quiescent 
SFR values. 

The SFRTOT luminosity function also differs signifi- 
cantly from that of the 2dFGRS luminosity function - a 
chi-squared test conducted as outlined above shows a van- 
ishingly small probability that the 2 luminosity functions 
are consistent, however, the SFRTOT sample differs in that 
it consists mainly of galaxies that are within ±1 mags of the 
value of Mj — 51og 10 h derived bv lNorberg et al.1 (|2002T l for 
the entire 2dFGRS sample. This can be understood by con- 
sidering the morphological analysis in section [3] where 69 
per cent of the sample were clas sified as spiral type galaxies 
and 29 per cent as S0/E type. Ijerien fc Tammannl (|l997I ) 
presented type-specific luminosity funcitons where the SO 
and spiral galaxy luminosity functions are best described by 
a Gaussian function. Thus, as expected, our Schechter func- 
tion fits have a > — 1, and the function behaves similarly to a 
modified Gaussian (see Table[2)|. The faint tail is supplied by 
a small number of low luminosity dwarf type galaxies where 
the luminosity function behaves like a Schechter function. 

It should also be noted that a burst of star formation 
will increase the 6j magnitude of a galaxy, hence pushing 
it towards brighter magnitudes up the luminosity function. 
The galaxies at the faint end of the SFRNORM sample 
would most likely fall below the magnitude cut off of the 
2dFGRS survey if they were not starbursting. 



5 LOCAL ENVIRONMENT 

We next examine the local environment of our starburst 
galaxies, and quantitatively assess whether it is of relevance 
to the starburst phenomenon. In doing so, we are particu- 
larly interested in the presence of near neighbour galaxies, 
the frequency of which was seen to be quite high when mor- 
phologically classifying our starburst samples visually (see 
Section 3). We use the SSS to extract catalogues of objects 
within 30 arcmin radius from our starburst galaxies. The 
approximate magnitude limits for the photographic plates 
scanned to produce the SSS are rp=21.5 and 6j=22.5. The 
objects in each catalogue are assigned a class parameter 1 
through 4 based on morphological parameters and areal pro- 
file. In this classification scheme, galaxies are class=l, stars 
class=2, unclassified class=3 and noise class=4. During the 
analysis, we use only those catalogue members with SSS 
class=l, noting there may be some contamination by stars. 

We conduct experiments designed to test for correla- 
tions between the starburst galaxy and its environment on 
scales less than 1 Mpc. During this analysis, we assume all 
neighbours lie at the same redshift as the starburst and con- 
vert angular separations to projected physical separations 
accordingly. Further to determining whether our starbursts 
have near neighbours, we would like to estimate the magni- 
tude of the gravitational perturbation due to the presence 
of the nearest neighbour. Hence we define an arbitrary mass 
ratio whereby if the near neighbour is more massive than 
one third of the mass of the starburst of interest, it con- 
stitutes evidence of some major interaction, whereas if the 
mass of the nearest neighbour is less than one third the mass 
of the starburst, it constitutes a minor interaction. By as- 
suming the mass-to-light ratio in the band is the same 
for the starburst and near neighbour, a mass ratio of 1:3 
constitutes a magnitude difference of 1.2 mag. We then de- 
fine bright and faint near neighbours (ie. major and minor 
interaction), by comparing the neighbour tf magnitude to 
that of the starburst magnitude corrected for the bright- 
ening due to the starburst, r F (see Appendix A for more 
details on how rp is derived), such that a faint neighbour 
has rJ±1.2<rF^ 21.5 and a bright near neighbour has 
r P +1.2. 

As an alternative local environment indicator, we use 
the surface density of the five nearest bright neighbours, 
E = 5/7rd| where ds is the transverse physical distance to 
the fifth nearest bright neighbour. 

The final test is designed to estimate the interaction 
strength due to the nearest neighbour. We use the Q- 
parameter, which estimates the ratio of the tidal (Fudai = 



GMnnRsbRs^) to binding force (Ft,, 
(|Daharilll984h hence 



GMsbR 



SB) 
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Figure 8. The luminosity distributions and the fitted Schechter luminosity functions for the SFRNORM (green), SFRTOT(orartge), 
SFRTOT+(pinfc) and 2dFGRS(6i«e) samples. The luminosity functions were derived using a stepwise maximum likelihood method and 
normalized to a constant source surface density of a = 124.6 deg -2 . The errors are fractional Poissonian errors. 



„ {Mnn\{D S b\ 3 (DnnDsb) 1 - 5 

Q(X {m^){r^) * IP ■ (6) 

Here, Msb / Mmn is the starburst/near neighbour mass and 
it is assumed that Mnn oc (Dnn) 1 ' 5 and Msb oc (Dsb) 1 ' 5 
jRubin et alJll982h , where L>atjv(= 2R nn ) is the nearest 
neighbour semi-major axis physical diameter, whilst Dsb(= 
2R$b) is the semi- major axis physical diameter of the star- 
burst of interest and S — R S e P where S is the projected phys- 
ical separation of the starburst and its nearest neighbour 
and R sep is the true three-dimensional separation. D$b and 
Dnn are derived from the intensity weighted semi-major 
axis given in the SSS catalogues which are used in conjunc- 
tion with the intensity weighted semi-minor axis to define 
an e lliptica l area for determination of isophotal magnitudes 
(see IStobid (1986) for more details). 

In order to gauge the significance of our results, for 
each test we compare our starburst sample results to that of 
a random sample containing 1000 galaxies. The selection of 
the random sample is crucial in order to negate systematic 
biases. Since we are using projected information to seek cor- 
relations with physical associations, it is important the ran- 
dom sample has the same redshift distribution as the star- 
burst sample of interest. This ensures both samples have the 
same systematic contributions due to foreground and back- 



ground galaxies. The random sample should also have the 
same magnitude distribution in each redshift bin. This is 
because we are using the rp magnitudes as a proxy for mass 
in an attempt to estimate the interaction strength. Hence 
we select separate random samples for the SFRNORM and 
SFRTOT samples, each with redshift distributions deter- 
mined by the starburst sample. The magnitude distributions 
were chosen to correspond to the brightening corrected mag- 
nitude starburst distributions. The approximate rp magni- 
tude limit for the 2dFGRS survey is 19 (assuming on 
average fcj-rF=0.5), hence the random sample contains only 
galaxies with tf ^ 19, whilst we only use starburst galaxies 
with brightening corrected magnitudes r F ^ 19. This addi- 
tional constraint on the starburst magnitude produced sam- 
ples of 197, 412 and 239 galaxies for the SFRNORM, SFR- 
TOT and SFRTOT+ samples, respectively. 

Figure [5] shows the near neighbour results for the 
SFRNORM sample. The errors presented are Poisson y~N 
errors. The Kolmogorov-Smirnov (KS) test was used to de- 
termine the probability that the random and starburst sam- 
ples were drawn from the same parent distribution. The KS 
statistic measures the maximum value of the absolute differ- 
ence in the cumulative probability distribution of two sam- 
ples and determines the significance of the difference from 
the KS-statistic distribution. Hence, we determine a prob- 
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Figure 9. Near neighbour analysis results for the SFRNORM sample. The upper left plot shows the distribution of distances to the 
nearest faint neighbour, the upper right plot shows the distribution of distances to the nearest bright neighbour, the lower left plot 
shows the distribution of the surface density of the nearest five bright neighbours, and the lower right plot shows the distribution of the 
logarithm of the Q-parameters as described in Section [5] The KS probabilities (see text) are 0.16, 0.015, 0.13 and 0.0017 for the faint 
near neighbour, bright near neighbour, surface density and Q-parameter, respectively. 



ability, p, that the KS statistic could exceed the observed 
value due to some random fluctuation in the data sets. A 
small value for p indicates two samples are not drawn from 
the same parent population. The p values determined for 
the faint nearest neighbour, bright nearest neighbour, sur- 
face density of bright neighbours and the Q-parameter are 
0.16, 0.015, 0.13 and 0.0017 respectively. These probabili- 
ties imply the bright near neighbour and Q-parameter dis- 
tributions differ significantly (at the 98.5 and 99.8 per cent 
levels, respectively) from the random distributions, with the 
bright near neighbour distribution having an excess of neigh- 
bours within ~ 20 kpc, whilst the Q-parameter distribution 
has a higher proportion of galaxies with log(Q) > — 1. We 
find 38 per cent of our SFRNORM galaxies have a bright 
neighbour within 20 kpc, and upon cross correlating with the 
morphological analysis in Section[3]an additional 16 per cent 
of this sample are classified as merger/interacting or possi- 
ble merger/interacting yet have no bright neighbour within 
20 kpc. 

Figure fTUl shows the near neighbour results for the SFR- 
TOT and SFRTOT+ samples, with Poissonian errors pre- 
sented. We used the KS test to derive p values for the faint 



nearest neighbour, bright nearest neighbour, surface den- 
sity of bright neighbours and the Q-parameter of T(.76), 
4.6 x 10" 5 (5.2 x 10" 5 ), 0.3(0.01) and 5 x 10~ 8 (6.7 x 10" 8 ), 
respectively, where the results for the SFRTOT+ sample are 
shown in brackets. Again, we find the bright near neighbour 
and Q-parameter distributions differ significantly from the 
random distributions (at > 99 per cent level), where the 
bright near neighbour distribution shows an excess of near 
neighbours within ~ 20 kpc and the Q-parameter distribu- 
tion has a higher proportion of galaxies with log(Q) > — 1. 
In addition to the 15 per cent for which we see a bright near 
neighbour within 20 kpc, we find 25 per cent are classified as 
merger/interacting or possible merger/interacting yet have 
no bright near neighbour within 20 kpc. Of the 22 per cent 
of the SFRTOT sample classified as E/S0, we find only one- 
quarter harbour a near neighbour within 50 kpc, hence we 
have a significant population of isolated, starbursting E/S0 
galaxies. 

The fact that the E5 distributions for the starbursts 
do not differ significantly f rom the overal random samples 
(similar to the findings of iMateus fc Sodrd (|2004l )). along 
with the statistically significant excess of bright near neigh- 
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Figure 10. Same analyses as shown in Fig.[9]for the SFRTOT (orange) and SFRTOT+ (blue) samples. The KS probabilities (see text) 
are 0.1 (0.76), 4.6 X 10~ 5 (5.2 X 10~ 5 ), 0.3 (0.01) and 5.0 X 10~ 8 (6.7 X 10~ 5 ) for the faint near neighbour, bright near neighbour, surface 
density and Q-parameter, respectively where the SFRTOT+ values are given in brackets. 



bours within 20 kpc for both starburst samples, suggest that 
a near neighbour plays an important role in triggering a 
starburst. These near neighbour results are consistent with 
previous results showing a neighbour within 50 kpc corre- 
lates with enhanced star formation rates (Nikolic et al.|[20o3 ; 
lLambas et al.ll2003l ; [Barton. Geller fc Kenvonl2000l ). The re- 
sults are also consistent with those of Freedman et al. (2006) 
who find that only galaxy pairs with small luminosity con- 
trasts show any trend of inc reased Ha emission at smaller 
projected radii. iGotol (|2005l ) also find a higher number of 
near neighbours within 50 kpc for 266 E+A galaxies when 
compared to field galaxies, consistent with our results. The 
Q parameter confirms the importance of a near neighbour, 
with both starburst samples showing excess galaxies with 
log(Q) > —1 compared to the random samples. Combining 
the near neighbour results with the morphological analysis 
in Section [3] we conclude that galaxy-galaxy interactions are 
important in triggering 40-50 per cent of the starbursts in 
our samples. Further to this, we note that there is no excess 
of faint near neighbours in either sample, implying that the 
interaction required to trigger a massive starburst needs to 
be significant. 



6 LARGE SCALE ENVIRONMENT 

Here we investigate the environments and clustering prop- 
erties of our starburst galaxy samples on more global scales, 
that is those larger than 1 Mpc. 



6.1 Membership of groups from the 2PIGG 
catalogue 

The 2PIGG galaxy groups catalogue (|Eke et al.ll2004al ) was 
compiled from the 2dFGRS catalogue using a friend-of- 
friends percolation algorithm. The catalogue was drawn 
from the contiguous SGP and NGP regions. For the pur- 
poses of our analysis, we restrict our attention to groups at 
redshifts z < 0.12, since above this limit interloper con- 
tamination increases, whilst the fraction of the observed 
group luminosity drops below half the total group lumi- 
nosity. There are two catalogues each for the NGP and 
SGP, the first contains the 2dFGRS galaxy information 
along with a flag of for ungrouped or the group ID 
number if the galaxy is linked to a group, the second 
contains information pertaining to group properties such 
as velocity dispersion, the unweighted number of group 
members and the rms projected galaxy separation. The 
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catalogues are publicly available and can be obtained at 
http://www.mso.anu.edu.au/2dFGRS/Public/2PiGG71 



We would like to ascertain whether our starburst galax- 
ies reside in particular types of groups when compared to 
a sample of 5000 random galaxies selected from the par- 
ent 2dFGRS catalogue such that their redshift distribution 
is the same as that of the starburst sample. There are a 
number of methods to determine group properties from the 
2PIGG catalogues. These include dynamical mass estimates, 
the numbe r of galaxi e s withi n each group and total group 
luminosity. lEke et all |2004bl ) showed that the weighted to- 
tal group luminosity is the most robust estimator of group 
size, since the dynamical estimators have large errors for 
the small groups, whilst the number of observed members 
residing in a higher redshift group can be significantly lower 
than the actual number due to the survey magnitude limit. 
Here we have used t he weighted obser ved group bj band 
luminosity defined bv lEke et al.l {2004b) as 



L b 



WiL 



i,bj , 



(7) 



total group lumi nosities by no more than ~ 10 per cent 



where Wi is the weight correcting for the incompleteness of 
the 2dFGRS, £;,6, is the bj luminosity of the ith group 
member and n ga i is the number of galaxies assigned to 
the group. From the observed weighted group luminos- 
ity, the total weighted group luminosity was determined 
by dividing by the incomplete Gamma function F(a + 
2, Lmi„/L*)/(r(a + 2)). This allows extrapolation from the 
lowest luminosity galaxy detectable at the group redshift, 
L m in, to zero luminosity, to correct for group members ly- 
ing below the survey magnitude limit. The correction is de- 
rived from the integral of the Schechter luminosity function 
in the luminosity r ange to L m i n , u sing global Schechter 
parameters given bv lEke et all (|2004bh as (M^ ,«)=(- 19.73, 
-1.18). The luminosity function is known to vary amongst 
groups of different masses, however using luminosity func- 
tion parameters derived on a group mass basis changes the 
up 

|Eke et al.ll2004bl ). 

In Figure [11] we plot the distributions of the corrected 
total group luminosities for the starburst and random sam- 
ples. We find 52, 58 and 56 per cent of the galaxies in 
the SFRNORM, SFRTOT, and SFRTOT+ samples reside 
in groups, respectively. Notably, these percentages are no 
different to those for the 2dFGRS galaxy population as a 
whole, with both random samples having 56 per cent of their 
galaxies residing in groups. Both the chi-squared and KS 
tests were used to determine the probability that the star- 
burst and random samples were drawn from the same par- 
ent distribution with both tests giving probabilities of less 
than 1 per cent of this occurring for each starburst sample. 
The distribution of luminosities for the SFRNORM sample 
are skewed to fainter and presumably poorer groups com- 
pared to the random sample, whilst the SFRTOT sample is 
strongly peaked at around 8L*, ie. medium sized groups. 



6.2 Membership of rich clusters 

|Pe Propris et all ()2002h have compiled a list of galaxy clus- 
ters from the Abell, APM and EDCC catalogues lying within 
the 2dFGRS regions. Redshifts, velocity dispersions and 



cluster centroids were measured for the majority of the clus- 
ters using the 2dFGRS data. We use both redshift informa- 
tion and cluster-centric distance to determine cluster mem- 
bership for our starburst galaxies, using only those clusters 
with measured redshifts. The velocity dispersion was used 
to determine the radius, T200, at which the density of the 
cluster exceeds that of the critical density by a factor of 200 
on the assumption the cluster is an isothermal sphere (see 
Carlberg et al. 1997). r2oo is an approximation of the virial 
radius of the cluster and is defined by: 

r200 = iWM' (8) 

where a v is the velocity dispersion, z is the cluster redshift 
and H(z) = Ho^Q m (l + z) 3 + Oa- For those clusters with- 
out a measured a v we used the average a v for the cluster 
sample in determining r2oo- We determine the peculiar ve- 
locity of the starburst with respect to the cluster, 



c((l + z pe 



1) 



((l + z pec ) 2 + l) 



(9) 



where z pec = (z s b — z)/(l + z) is the non-cosmological red- 
shift contribution of the starburst galaxy to its observed 
redshift, z 3 b, due to its peculiar velocity within the clus- 
ter. We determined r, the projected physical distance of the 
starburst from the cluster centroid as r = A9x{z)/(1 + z) 
where x( z ) is the radial comoving distance to the cluster 
and A8 is the angular separation of the cluster centroid and 
the starburst galaxy. Those starburst galaxies with r ^ T2oo 
and Vpec ^ 3a v are classed as cluster members, whilst those 
with r2oo < t ^ 2r2oo and v pec ^ 3a v are classed as lying in 
the cluster outskirts and those with 2r2oo < r ^ 4r2oo and 
v pec ^ 3cr„ are classed as being within the infall regions of 
the cluster. 

Using the above definitions, we find for the SFRNORM 
sample that 3, 10 and 19 per cent of starbursts lie within the 
cluster, outskirt and infall regions, respectively, compared 
to 8, 8 and 21 per cent respectively, for our random sam- 
ple of 5,000 galaxies. The same analysis performed on the 
SFRTOT (SFRTOT+) sample produces 5 (5), 6 (7) and 20 
(22) per cent within the cluster, outskirt and infall regions, 
respectively, compared to 9, 10 and 24 per cent for our ran- 
dom sample of 5,000 galaxies. Hence, we conclude that the 
starburst galaxies within our sample are not preferentially 
found within clusters, with ~ 30 per cent total found within 
4 x r2oo Mpc of a cluster centroid. 



6.3 Residence within large-scale overdensities 

As another method of quantifying the large scale structure 
surrounding our starburst galaxies, we measured the over- 
densities, S, within comoving spheres of radii 1 ^ ro ^ 
15 Mpc in 1 Mpc intervals. Specifically, we used the comov- 
ing distance of each surrounding galaxy to transform its 
right ascension, declination and redshift into 3-D cartesian 
coordinates, making no attempt to correct for peculiar ve- 
locity distortions in redshift space, and measure its physical 
separation from the starburst galaxy of interest. For each 
comoving sphere centred on a starburst galaxy, we add up 
the number of galaxies in the 2dFGRS catalogue within ro 
giving the number of observed galaxies within the sphere, 
N bs. The number of galaxies expected in each comoving 
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Figure 11. The distribution of corrected total group luminosities for the SFRNORM sample {left) and SFRTOT sample {right) compared 
to the same distributions for the associated random samples (see text for details). Error bars show Poissonian yN values, where N is 
the number of galaxies in each bin. Both starburst samples differ from their respective random samples at greater than 99 per cent 
significance. 



sphere in the absence of clustering, N exp , was determined 
using random mock catalogues generated by the publically 
available 2dFGRS mask software written by Peder Norberg 
and Shaun Cole (see |http://www2. aao.gov.au/2 dFGRS/ ). 
Again, using 3-D cartesian coordinates, we measure the 
physical separation between the mock catalogue galaxy and 
the starburst of interest. We obtain N exp by adding up the 
number of mock galaxies surrounding each starburst galaxy 
within Vq. The average overdensity within each comoving 
sphere is calculated as follows: 



(V30 



i=l 



Not 
N ex 



(10) 



where n s i is the number of starburst galaxies used in the 
analysis. The random mock catalogue was generated such 
that the sample size is 30 times that of the 2dFGRS cata- 
logue out to z=0.18 with a limiting magnitude of 6j=19.45, 
hence the multiplicative factor of 30 in eqn. 1101 The larger 
sample size for the mock random catalogue was used in or- 
der to minimize the noise in the denominator caused by low 
number statistics. 

We compared to a random sample selected from the 
2dFGRS to have the same redshift distribution as the star- 
burst sample of interest, such that redshift dependent biases 
are minimized. We used only the SGP and NGP regions for 
this analysis, the results of which are presented in Figure [121 
The error bars presented are the variance, a/^fn^, of the 
overdensity distribution where a is the standard deviation 
of the distribution of mean overdensities and n s b is the num- 
ber of starburst galaxies in the sample. The results clearly 
indicate that both starburst samples reside in lower density 
environments than the randomly selected 2dFGRS galax- 
ies on all scales from 1 Mpc to 15 Mpc. A common measure 
of environment is the overdensity within spheres of radius 
8/i -1 Mpc. On this scale, we measure average overdensities 
of 0.92 (±0.08) and 0.84(±0.06) for the SFRNORM and 
SFRTOT samples, respectively, which are both consistent 
(within errors) with that measured by iBlake et al.l (|2004l ) 
for the 2dFGRS E+A sample. For the random samples we 



measure 1.17(±0.02) and 1.19(±0.02) for the SFRNORM 
and SFRTOT random samples, respectively. Despite using 
a different method to measure the number of galaxies ex- 
pected in each sphere in an unclustered Universe, the results 
are, within the errors ? consistent with th e value measured for 
a random sample by IBlake et al.l (|2004l ) . 



6.4 Spatial cross-correlation function 

The spatial cross-correlation function measures the cluster- 
ing properties of galaxies by estimating whether a particular 
type of galaxy is biased to inhabit high mass or low mass 
haloes when compared to the bias of another type of galaxy. 
If we approximate a linear bias for the galaxies in the star- 
burst samples, b s b, then their spatial auto-correlation func- 
tion, £ s b,sb, as a function of spatial separation, r, takes the 
form 



£ s b,sb{r) = & si> £ m , m (r), 



(11) 



where £ m , m (r) is the spatial auto-correlation function of the 
underlying mass density field. 

Due to the relatively low number of galaxies in our star- 
burst samples, shot noise dominates the auto-correlation 
function, hence we measure the cross-correlation function, 
£sb, g , of the starburst galaxies with the rest of the 2dFGRS 
catalogue. In the linear bias approximation, 



£sb, g (r) = & s fcb g ^ m , m (r), 



(12) 



where b g is the bias factor for an average 2dFGRS galaxy. 
We can estimate the relative bias of the starbursts to the 
2dFGRS galaxies by measuring the auto-correlation func- 
tion for the 2dFGRS galaxies, £9,3 ( r ) = b^ m .m{r), such 
that 



bsb 

ba 



£sb,g(r) 



(13) 



We have measured all correlation functions in redshift 
space, making no attempt to correct for peculiar velocities. 
The cross-correlation function is measured by comparing 
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Figure 12. The average overdensity surrounding each starburst galaxy within comoving spheres of radius ro obtained by comparing the 
number of galaxies in the 2dFGRS within ro with that of a mock random catalogue. The analysis was repeated on a sample of 10,000 
random 2dFGRS galaxies selected such that the redshift distribution is the same as the starburst sample of interest. The figure on the 
left shows the results for the SFRNORM sample, whilst that on the right shows the results for the SFRTOT and SFRTOT+ samples. The 
error bars represent the variance in the overdensity distribution; for the random samples these are negligible and hence not presented. 



the cross-pair counts of the starburst samples with the full 
2dFGRS catalogue and an unclustered random distribution 
with the same redshift distribution and selection mask as 
the 2dFGRS sample, whilst in order to minimize fluctua- 
tions due the selection function of the starburst galaxies, we 
measure the cross pair counts of the 2dFGRS with an un- 
clustered random sample generated with the same redshift 
distribution as the starburst sample of interest. The ran- 
dom distributions were generated using the same publicly 
available 2dFGRS mask software as in Section T6. 31 and con- 
tain the same number of galaxies as the real data samples 
they ar e simulating. We use a modified Landy-Szalay esti- 
mator l|Blake et alj|2006h for the cross- correlation function, 
£sb,g(s), which is given by 



N ab ,g{s) N sblR g(s) Ng, Rsb (s) 



NRsb,Rg(s) N Rsb:Rg (s) N Rab , Rg (s) 



+ 1, (14) 



where N sbig (s) is the starburst-2dFGRS pair cross count, 
N sbtRg (s) is the starburst-simulated unclustered random 
2dFGRS sample pair cross count, N gtRsb (s) is the 2dF- 
GRS -simulated unclustered random starburst sample pair 
cross count, N RsbtRg is the cross pair count for the two 
random samples and the s denotes a redshift space sepa- 
ration. Each of the N sb , R g(s), Ng, Rsb (s) and N Rsb , Rg are 
measured by generating 10 random catalogues and tak- 
ing the average of the 10 measurements. In order to esti- 
mate the auto-correlation function for the 2dFGRS galaxies, 
£q, q , for comparison, w e use the Landy & Szalay estimator 
i Landv fc Szalavlll993h . given by 



£ 9 , 9 ( s ) = 



N g ,g(s) Ng, Rg ( S ) 



N R g, R g(s) N R g, R g(s) 



(15) 



where N BtB is the auto-pair count for the 2dFGRS 
galaxies, N g , Rg is the 2dFGRS -random pair cross count and 
N RgtRg (s) is the random-random pair auto count. 

We include the contiguous NGP and SGP strips in our 
analysis, but not the random fields since the variance in the 
correlation function is large due to edge effects. The results 
are displayed in Figure 1131 We do not assign Poisson error 



bars here, since they are known to underestimate the true 
variance of the estimators in equation s ll4l and ll5l bv a signif- 
icant amount (|Landv fc Szalavlll993l ). The errors presented 
are estimated using the 'jack-knife' approach whereby the 
NGP and SGP strips are divided into forty regions and the 
correlation function estimation is repeated forty times, each 
time leaving out one region. The error for each separation 
bin is then estimated by multiplying the resulting standard 
deviation across the forty subsamples by y/40. 

The results are presented in Figure [13] where £ s , 9 (s) 
and £ sb ,g(s) are plotted for the SFRNORM, SFRTOT, SFR- 
TOT+ and the 2dFGRS catalogues. We also plot the ra- 
tio of U,s(s)/£s,s( s )- The SFRNORM, SFRTOT and SFR- 
TOT+ samples are much less clustered than the entire 2dF- 
GRS sample, on all scales (apart from the smallest separa- 
tion bin for the SFRTOT and SFRTOT+ samples, which are 
consistent with the 2dFGRS sample within errors) imply- 
ing that b sb < b g . In fact, the average ratio of b ab /b g is 0.53, 
0.69 and 0.73 for the SFRNORM, SFRTOT and SFRTOT+ 
samples, respectiv ely. These results are consistent with the 
Blak e et al.l (|2004l ) ' average Balmer' E+A catalogue, which 
is found to be somewhat less clustered than the 2dFGRS 
catalogue, although it is noted their sample size is only 50, 
hence the results are tentative. The resu lts are also consis- 
tent with those o f iMadewick et all ||2003D . who find that ac- 
tive galaxies, defined by the r\ parameter (jMadgwick et all 
2002), are less clustered than passive galaxies out to real 
space separations of 10/i _1 Mpc. 



7 SUMMARY AND DISCUSSION 

We have selected two samples of starburst galaxies, each con- 
taining 418 galaxies, using different methods for deriving the 
star formation rate from the Ha emission line. The selec- 
tion methods are complementary, in that selecting starburst 
galaxies based on W R a alone is known to be biased towards 
low mass galaxies (our SFRNORM sample), whereas the se- 
lection of galaxies by measuring the total absolute star for- 
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Figure 13. Spatial correlation function measurements for the 
SFRNORM, SFRTOT, SFRTOT+ and 2dFGRS samples. The 
plotted errors are derived using a jack-knife re-sampling tech- 
nique, as described in the text. The lower plot shows the ratio 
of Lb, g (s)/tg,g(s) for the SFRNORM, SFRTOT and SFRTOT+ 
samples. 

mation rate selects a sample of more luminous, and presum- 
ably more massive galaxies (our SFRTOT sample). Hence, 
we cover a broad spectrum in mass with our two starburst 
samples. We study the environments and morphologies of 
the galaxies in these samples in an effort to ascertain mech- 
anisms which may be important in triggering a burst of star 
formation, and in truncating the star formation leading to 
the E+A type spectrum. 

7.1 Summary of results 

The morphological analysis of our starburst galaxies indi- 
cates: 

(i) Both starburst samples contain galaxies covering the 
full range of Hubble types, with early types (E/SO, SE) being 
more dominant in the SFRTOT sample than the SFRNORM 
sample. 

(ii) 20-30 per cent of the starburst galaxies, irre- 
spective of how they are selected, exhibit signs of major 
merger/interaction or tidal activity. 

(iii) 29 per cent of the SFRTOT sample are classified as 
E/SO and of these early types 69 per cent show no evidence 
for merging/interaction nor do they harbour a bright near 
neighbour within 20 kpc. 

We have conducted the same environmental analyses 
on a sample of random 2dFGRS galaxies, selected such that 
systematic effects relevant to the particular analysis are min- 
imized. Comparison of our starburst samples with these ran- 
dom samples shows: 

(i) The luminosity function of the SFRNORM sample 
has a significantly steeper faint-end (Mb > — 18) slope com- 
pared to the 2dFGRS catalogue. As noted above, selection 
based on Who, alone is biased towards selecting low mass 
galaxies, hence we should expect a high number of low lu- 
minosity galaxies in our sample. In contrast, the SFRTOT 
sample shows a dearth of galaxies at M& > —18, and is dom- 
inated by L* and brighter galaxies. The dearth of fainter 
galaxies is presumably caused by the fact that the low lumi- 



nosity galaxies don't have high absolute SFRs, but do have 
high SFRs compared to their quiescent states. Indeed, the 
majority of galaxies in the SFRNORM sample have SFR- 
TOT values < 10 M Q yr~ x . 

(ii) Both samples have a significant overdensity of bright 
near neighbours within 20 kpc. The SFRTOT sample also 
shows an excess number of galaxies experiencing tidal forces 
due to a near neighbour. This excess is also present in the 
SFRNORM sample, although at a lower level of significance. 
For the local galaxy surface density and projected distance 
to faint near neighbour studies, we find the distributions for 
the starburst samples are consistent with that of a random 
sample. 

(iii) On larger scales, we find those SFRNORM galax- 
ies residing in groups preferentially inhabit low luminosity 
groups, whilst a higher fraction of those SFRTOT galaxies 
residing in groups inhabit average luminosity groups. We 
find our starbursts do not preferentially inhabit clusters of 
galaxies. Upon measuring the overdensity within comoving 
spheres on scales 1-15 h~ Mpc we find that starbursts lie 
in underdense regions compared to random samples on all 
scales. The cross-correlation function also shows that our 
starburst galaxies are less clustered on scales larger than 
2h~ 1 Mpc compared to the 2dFGRS galaxies. 

(iv) We find throughout the analysis that the SFR- 
TOT+ sample shows no significant differences from that of 
the SFRTOT sample. 

From this analysis, we can conclude that starburst 
galaxies are less clustered and hence less biased on large 
scales when compared to passive galaxies and appear to be 
mainly found within the field rather than clusters of galax- 
ies. Local environment appears to be the key determinant in 
triggering a starburst, in the form of mergers and tidal inter- 
actions. The analysis of the types of groups our starbursts 
inhabit is consistent with the merger/interaction hypothesis 
when we consider that the SFRNORM sample is made up of 
primarily low luminosity /mass galaxies with presumably low 
velocity dispersions, whilst the SFRTOT sample consists of 
more luminous/massive galaxies which should have higher 
velocity dispersions. Mergers are rare in environments where 
the galaxy velocity dispersion is much higher than the stel- 
lar velocity dispersion, hence if mergers are the dominant 
mechanism triggering a starburst then we would expect to 
find low mass starbursts in low mass groups, and high mass 
starbursts in intermediate mass groups, consistent with the 
results presented here. 

It is also clear from our analysis that a large fraction 
of our starbursts show no evidence for merger/interaction, 
nor do they have a bright close near neighbour. Hence, some 
form of galaxy interaction can not have triggered the star- 
burst. What mechanism is triggering the starburst in these 
galaxies? Also of note is that 22 per cent of the SFRTOT 
sample is classified as E/SO with no signs of merger activ- 
ity, whilst our near neighbour analysis shows only 25 per 
cent of these undisturbed E/SO galaxies harbour a bright 
near neighbour within 50 kpc. It may be that there is some 
merger activity, but either our criterion is too strict to de- 
fine the galaxy as a merging/interacting system, or the SSS 
imaging is simply too poor to resolve any evidence of this. It 
may also be that there has been some tidal interaction with 
a near neighbour some time ago, and the near neighbour 
is now at some distance from the starburst, hence does not 
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show up in our near neighbour analyses. For example, if a 
starburst was triggered by a galaxy interaction 500 Myrs ago 
and the near neighbour galaxy was traveling at 100 km s -1 , 
then the near neighbour would be some 50 kpc distant from 
the starbursting galaxy. This oversimplified example serves 
to show how evidence for a merger or tidal interaction may 
be missed in our analyses. There is also the possibility that 
these galaxies are in the later stages of swallowing a gas-rich 
faint dwarf galaxy which is no longer distinguishable from 
the host galaxy. However, if these starbursts were not trig- 
gered by some interaction, then an internal mechanism may 
be required to explain the burst. Detailed follow-up spatially 
resolved spectroscopy and high resolution imaging is critical 
in order to determine the nature of the triggering mechanism 
in these galaxies. 

7.2 Comparison with 2dFGRS E+A galaxies 

Concerning starburst galaxies as E+A progenitors, it is 
important to ask if our samples are likely to evolve into 
E+A galaxies bas ed on the evidence presented here and in 
Blake et al.l d2004h. It is important to note at the outset here, 
that iBlake et all (jioolh found only 56 'true' E+A galaxies, 
whilst each of our starburst samples contain 418 galaxies de- 
spite probing a smaller volume. Perhaps only a small fraction 
of our starburst galaxies will go through an E+A phase. 



7.2.1 Luminosity functions 

Comparison of the E+A and SFRNORM luminosity func- 
tions shows a marked difference at the faint end (M;,j — 
51og 10 h, < —18) from which it might be inferred that our 
SFRNORM sample galaxies will not evolve into E+A galax- 
ies. However, if our sample is predominantly dwarf star- 
bursts which undergo periodic bursts of star formation with 
a duty cycle of 10 Myr, then we would expect post-starburst 
dwarf galaxies would display the characteristic E+A spec- 
trum on very short time scales. Coupled with the fact that 
a large fraction of our SFRNORM sample are close to the 
magnitude limit of the 2dFGRS survey and would not be 
detected in their un-brightened, quiescent state, we expect 
low luminosity E+A galaxies would be v ery rare in the 2dF - 
GRS survey. Indeed, this is seen in the IBlake et all (|2004l ) 
survey. Thus, it is unlikely that ou r SFRNORM samp le will 
evolve into the E+As observed by IBlake et all (|2004f ). 

Co mparison our SF RTOT luminosity function with 
that of IBlake et all (12004 ) ' average Balmer' sample (EW ~ 
4.5 A), shows that the two samples differ at 97.7 per cent 
significance level in the range —17.75 > Mb, — 51og 10 h, > 
—20.75 where the difference occurs due to a higher num- 
ber of M bj - 51og 10 /i = -19 ± 1 galaxies in the SFRTOT 
sample. Taking this at face value, it would seem that the 
our starburst galaxies can not be progenitors to E+A galax- 
ies. However, since we are measuring the bj band luminos- 
ity function, we would expect the luminosity function of a 
starburst sample to be significantly different from that of a 
sample of quiescent galaxies with the same mass distribu- 
tion. This is because the radiation emitted by the starburst 
is dominated by a population of young, massive O-type stars 
which emit predominantly at blue wavelengths, and hence 
the bj luminosity of a starburst galaxy is significantly en- 
hanced compared to its quiescent state. After the starburst 



is truncated and the galaxy enters the E+A phase, its bj 
luminosity gradually fades back towards its quiescent value. 
This evolution in bj magnitude means the luminosity func- 
tion is not a particularly good discriminator for the two 
samples. The k-band luminosity function would provide a 
much better discriminator here, since it is less affected by 
the light emitted from starburst regions. Therefore, the dif- 
ference in the bj luminosity functions alone is not enough 
to rule out an evolut ionary link betwee n our starbursts and 
the E+A galaxies of IBlake et"ai] |2004h . 



7.2.2 Local environment 



IBlake et al. I (|2004l ) find no statistically significant differences 
between their E+A samples and random 2dFGRS samples 
in their near neighbours analysis, whereas here we find clear 
differences, particularly for the bright nearest neighbours 
within ~ 20 kpc. This is also somewhat expected, since sim- 
ulations have shown that the E+A phase generally occurs at 
a late stage in mergers when t he cores of the two g alaxies can 
no lo n ger b e distinguished (B ekki et al.l bOO lb'), although 
iGotd lj2005h find an enhancement in the number of near 
neighbours in their E+A sample. Also, IBlake et alj (|2004h 
find a significant portion of their E+A sample exhibit tidal 
tails and merger morphologies and the sample is made up of 
predominantly E/S0 galaxies, hence they conclude that ma- 
jor mergers are an important formation process for E+As. 
We then conclude that a sub-sample of our starbursts are 
likely to evolve into E+A galaxies after the rapid star for- 
mation has consumed the available gas. 



7.3 Caveats 

It is worth noting some caveats to our conclusions. The 
first is the 2dFGRS fibre size, which is 2arcsec in diame- 
ter. This corresponds to 4kpc at z=0.11, hence the covering 
fraction of the fibre is in general much smaller than the 
average galaxy size. Studies of aperture effects have shown 
that measured properties such as corrected SFR and metal- 
licity are affected when the covering fraction of a 2arcsec 
fibre is < 20 per cent, which corresponds to a redshift of 
z=0.06 for galaxies with si milar properties to t hose in the 
Near Field Galaxy Survey (jKewlev et al.l l2005h . Thus we 
should expect our SFRTOT sample to have SFR overesti- 
mated for z < 0.06. This effect is somewhat nullified by 
the fact that we select our sample from the most extreme 
galaxies per redshift bin. The effect should also be negligi- 
ble in the SFRNORM sample due to two main reasons. The 
first being this sample is dominated by dwarfs, and indeed 
the distribution of diameters derived from SSS parameters 
peaks at ~ 3 kpc, hence the covering fraction drops below 
20 per cent at around z=0.035, affecting some 35 per cent 
of our SFRNORM sample. The second is our SFRNORM 
measurements are based on Whch, which is a relative quan- 
tity for which we see minimal trends out to z = 0.18. The 
second caveat comes from the fibre positioning in the 2dF- 
GRS survey. The fibres are most likely to be placed at the 
approximate centre of the galaxy, thus if the covering frac- 
tion is small, we only receive light from the central regions 
meaning we are biased to selecting nuclear starburst galax- 
ies. This may mean we miss particular types of starburst 
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galaxies which have non-nucleated starbursts, eg., those with 
bursts triggered by ram-pressure from the ICM where the 
burst is expected in the galaxy outskirts, or those where the 
burst is triggered by gas sloshing. Follow up long-slit spec- 
troscopy of a sub-sample of starburst galaxies will allow us 
to decipher how much these aperture affects should bias our 
results. 

7.4 Outlook 

Clearly, from this study it can be concluded that a sig- 
nificant fraction of starbursts are triggered by some form 
of merger/interaction event. However, equally important 
are those galaxies that show no evidence for interactions, 
and hence what might have triggered their starburst. Fol- 
low up studies of these galaxies including high resolution 
optical imaging (for morphological analysis), spatially re- 
solved spectroscopy (for mapping out the star forming re- 
gions across the galaxy and measuring stellar kinematics) 
and HI imaging (for constraining gas dynamics) will allow 
us to determine whether the triggering mechanism is inter- 
nal or external in these starbursts. 
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abundances of 20 and 60 per cent of the solar value for the 
SFRNORM and SFRTOT samples, respectively. Hence we 
use the corresponding metallicities in the BC03 models. We 
find that, on average, a galaxy in the SFRNORM sample 
is brightened by ~ 0.83 mag and a galaxy in the SFRTOT 
sample is brightened by ~ 0.44 mag in the r-p band due to 
a starburst. Thus, for the SFRNORM sample r F =r F +0.83 
and for the SFRTOT sample, r P =r P +0.44. 



APPENDIX A: CORRECTING FOR THE 
BRIGHTENING IN THE r F MAGNITUDE DUE 
TO A STARBURST 

During the near neighbour analysis, we assumed a constant 
mass to light ratio for the starburst and near neighbour 
galaxies. This assumption is subject to errors caused by the 
brightening of the tf magnitude due to the burst of star 
formation occurring in the starburst, hence, we attempt to 
correct for this using the stellar population synthesis models 
of Bruzual & Chariot (2003). Briefly, we superpose a 300 Myr 
burst of constant star formation onto an underlying 11 Gyr 
old stellar population which has an exponentially decaying 
star formation rate with an e-folding time of r = 3 Gyrs be- 
ginning 1 Gyr after the Big Bang. There are two main inputs 
for the models which we can derive from our spectra - the 
burst strength (ie. the ratio of the mass of stars formed in 
the burst to the underlying stellar mass) and the metallicity. 

The burst strength was determined from the Who) of 
the starburst sample of interest (for the SFRNORM the 
median W Ha = 145 A and for the SFRTOT sample the 
median Wna = 67 A. The Wh c correlates strongly with 
the 'b' parameter (IScalol I1986T ) which measures the ratio 
of the current star formation rate to the mean past star 
formation rate (IKennicutt et alj|l994h ■ Using Table 1. from 
iKennicutt et al.l (|l994l ) , we estimate the equivalent b param- 
eters for the median W Ha are 2.5 and 1 for the SFRNORM 
and SFRTOT samples, respectively. Multiplying the 'b' pa- 
rameter by the ratio of the assumed ages of the starburst 
(300 Myr) and underlying stellar population (11 Gyr), we 
can derive the burst strength. We derive burst strengths of 
7 and 3 per cent for the SFRNORM and SFRTOT samples, 
respectively. 

To determine the metallicity, we obtain the gas phase 
oxygen abundances using the EWs of 0111(4959,5007), 
[Oil] [3737) and H/3 to derive the R23 parameter (see Kobul- 
nicky et al. 2003) which in turn is used to determine the oxy- 
gen abundance, 12+log(0/H) (see Kobulnicky et al. 1999). 
The degeneracy between the R23 parameter and its relation 
to the oxygen abundance is broken by defining metal rich 
galaxies as those with log([NII] A6584/Hq)^ —1 and metal 
poor as those with log([NII]A6584/Ha)< —1 (see Lamaraille 
et al. 2004 and references t herein for a more thorou gh discus- 
sion). Eqns. (8) and (9) of iKobulnickv et al. ( 1999) are then 
used to determine the 12+log(0/H). We find mean oxygen 



